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Abstract — The design of modulation schemes for the physical 
layer network-coded two-way MIMO relaying scenario is con- 
sidered, with n-R antennas at the relay R, iia and Ub antennas 
respectively at the end nodes A and B. We consider the denoise- 
and-forward (DNF) protocol which employs two phases: Multiple 
access (MA) phase and Broadcast (BC) phase. It is known for the 
network-coded SISO two-way relaying that adaptively changing 
the networking coding map used at the relay, also known as 
the denoising map, according to the channel conditions greatly 
reduces the impact of multiple access interference which occurs 
at the relay during the MA phase and all these network coding 
maps should satisfy a requirement called the exclusive law. The 
network coding maps which satisfy exclusive law can be viewed 
equivalently as Latin Rectangles. In this paper, it is shown that for 
MIMO two-way relaying, deep fade occurs at the relay when the 
row space of the channel fade coefficient matrix is a subspace 
of a finite number of vector subspaces of C" A+ " B which are 
referred to as the singular fade subspaces. It is shown that 
proper choice of network coding map can remove most of the 
singular fade subspaces, referred to as the removable singular 
fade subspaces. All these network coding maps are obtainable by 
the completion of partially filled Latin Rectangles. For 2 A -PSK 
signal set, the number of removable and non-removable singular 
fade subspaces are obtained analytically and it is shown that 
the number of non-removable singular fade subspaces is a small 
fraction of the total number of singular fade subspaces. The 
Latin Rectangles for the case when the end nodes use different 
number of antennas are shown to be obtainable from the Latin 
Squares for the case when they use the same number of antennas, 
irrespective of the value of jir. For 2 A -PSK signal set, the singular 
fade subspaces which are removed by the conventional XOR 
network code are identified. Also, using the notions of isotopic 
and transposed Latin Squares, the network coding maps which 
remove all the removable singular singular fade subspaces are 
shown to be obtainable from a small set of Latin Squares. 

I. Background and Preliminaries 

We consider the two-way wireless relaying scenario shown 
in Fig[T] with multiple antennas at the nodes, where data 
transfer takes place between the nodes A and B with the help 
of the relay R. It is assumed that all the three nodes operate 
in half-duplex mode, i.e., they cannot transmit and receive 
simultaneously in the same frequency band. We consider the 
denoise-and-forward (DNF) protocol originally introduced in 
JT), which consists of the following two phases: the multiple 
access (MA) phase, during which A and B simultaneously 
transmit to R and the broadcast (BC) phase during which R 
transmits to A and B. Network coding map, which is also 
referred to as the denoising map, is chosen at R in such a way 
that A (B) can decode the messages of B (A), given that A 
(B) knows its own messages. 
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(b) BC Phase 
Fig. 1. The Two Way Relay Channel 



The concept of physical layer network coding has attracted 
a lot of attention in recent times. The idea of physical layer 
network coding for the two way relay channel was first intro- 
duced in [2], where the multiple access interference occurring 
at the relay was exploited so that the communication between 
the end nodes can be done using a two phase protocol. Infor- 
mation theoretic studies for the physical layer network coding 
scenario were reported in J3), |4). A differential modulation 
scheme with analog network coding for two-way relaying was 
proposed in Q. The design principles governing the choice 
of modulation schemes to be used at the nodes for uncoded 
transmission were studied in [6 |. An extension for the case 
when the nodes use convolutional codes was done in Q. A 
multi-level coding scheme for the two-way relaying scenario 
was proposed in J8). Power allocation strategies and lattice 
based coding schemes for two-way relaying were proposed in 

m. 

It was observed in J6) that for uncoded transmission, the 
network coding map used at the relay needs to be changed 
adaptively according to the channel fade coefficient, in order 
to minimize the impact of the Multiple Access Interference 
(MAI). For the single antenna two-way relaying scenario, 



a computer search algorithm called the Closest-Neighbour 
Clustering (CNC) algorithm was proposed in |6l to obtain the 
adaptive network coding maps resulting in the best distance 
profile at R. An extension to MIMO two-way relaying using 
the CNC algorithm was made in ifTOl . With a single antenna 
at the nodes, the MAI becomes severe for channel fade states 
referred to as the singular fade states lfl2l . An alternative 
procedure to obtain the network coding maps, based on the 
removal of singular fade states using Latin Squares was 
proposed in Ifl2l . A quantization of the set of all possible 
channel realizations based on the network code used was 
obtained analytically in |fl3ll . 

In this paper, it is shown that for the MIMO two-way 
relaying scenario, the MAI becomes severe when the row 
space of the channel fade coefficient matrix is a subspace of 
a finite number of vector subspaces of C™ A+ ™ B referred to 
as the singular fade subspaces. The notion of singular fade 
subspaces subsumes in it as a special case the notion of 
singular fade states introduced in [ 12] for the single antenna 
two-way relaying scenario. It is shown that the network coding 
maps totally avoiding the removable singular fade subspaces 
can be obtained by the completion of partially filled Latin 
Rectangles. 

Notations: Let CAf(0, a 2 I n ) denote the circularly symmet- 
ric complex Gaussian random vector of length n with covari- 
ance matrix a 2 I n , where is the identity matrix of order n. 
The binary field is denoted by F 2 . All the vector spaces and 
vector subspaces considered in this paper are over the complex 
field C. Let span(c!, C2, . . . Cl) denote the vector space over 
C spanned by the complex vectors ci, C2, . . . cl G C". For a 
matrix A, A T denotes its transpose. For a vector subspace V of 
C", V 1 - denotes the vector subspace {x : x T v = OVw G V}. 
For a matrix H G C mx ", 71(H) denotes the row space of H. 
The all zero vector of length n is denoted by 0„. For a vector 
x of length n, Xi, 1 < i < n denotes the i th component of 
x. For vector subspaces V\ and V2 of C", V\ < V2 means 
that V\ is a subspace of V2 and V\ ^ V2 means that V\ 
is not a subspace of V%. By ua x n B system, we refer to 
the two-way relaying system with ua and ng antennas at 
the nodes A and B respectively, with no restriction placed 
on the number of antennas at the relay n R . For a vector x, 
x [i-.j] j 1 5; i < 3 5; ti denotes the vector obtained by taking 
only the i th to j th components of x. 

A. Signal Model 

Consider the MIMO two-way relaying system with ua-, n R 
and n B antennas at the nodes A, R and B respectively as 
shown in Fig. [TJ It is assumed that the complex numbers 
transmitted in each one of the antennas at the end nodes belong 
to the signal set S of size M = 2 A , where A is an integer. 
Assume that A (B) wants to transmit a An^-bit (Ans-bit) 
binary tuple to B (A). At node A (B), the Xua {Xns) bits 
are spatially multiplexed into ua (ns) streams with each one 
of the streams consisting of A bits. The A bits in each one of 
the streams are mapped on to the signal set S using the map 
li : —> S and are transmitted. 



Multiple Access (MA) phase: Throughout, we assume 
a block fading scenario with the Channel State In- 
formation (CSI) available only at the receivers. Let 

X A = [Ll(s Al ),n( s A 2 ),-,n{sA nA )} T G S nA , X B = 

[^(ssj, ll(sb 2 ), m( s s„ b )] T G S ns denote the com- 
plex vectors transmitted by A and B respectively, where 
s Al , s A2 , —, s A „ A , s Bl , Sb 2 , -, s Bng £ F£. The received sig- 
nal at R is given by, 

Y R = H A x A + H B x B + Z R , 

where Ha of size n R x Ua and Hb of size n R x are 
the channel matrices associated with the A-R and B-R links 
respectively. The additive noise vector Z R is assumed to be 
CAf(0,a 2 I 1lR ). 

Let Sr(Ha, Hb) C C nR denote the effective constellation 
seen at the relay during the MA phase, i.e., 

S R (H A , H b ) = {H A x A + H B x B \x A £ S ,lA , x B £ S riB } . 

Let d m i n (HA, Hb) denote the minimum distance between 
the points in the effective constellation S r (Ha, Hb), i.e., 

d m in{H A ,H B )= min H A (xa — x' A ) + H B (x B - x' B ) \\ 

(x A ,x B ),(x' A ,x' B )eS n A + "B 

{ x a^b)^( x 'a' x 'b) 

(1) 

From ([T), it is clear that there exists values of Ha,Hb for 
which d mm (H A , H B ) = 0. 
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Fig. 2. Diagram showing the effective constellation at the relay for Hr f» 
[1+j 2 3 i + j] T 

Example 1: Consider the case when = ng = 2 
and n R = 1 and BPSK signal set {±1} is used at the 
end nodes. When [H A H B ] ~ [1 + j 2 3 4 + j] T , the 
effective constellation S R is as shown in Fig. [2] It is clear 
from Fig. E] that for [H A H B ] = [1 + j 2 3 4 + j] T , 

dmin(HA, H B ) =0. 

Let AS denote the difference constellation of the signal set 
S, i.e., AS = {s i -s' i \si,s' i €S}. 

Let us define Ax = \{xa ~ x' A ) T (x B — x' B ) T ] T G 
AS nA+nB , where Ax A = x A - x' A G AS nA and Ax B = 
x B - x' B G AS" b , and H R = [H A H B }. Also let H R = 
[hi hi ...h nii ] T , where the row vector h^,l < k < n R of 



length riA + n B is the fc th row of Hr. From (Q]), it follows 
that, 
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From (O it follows that d. 
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that fc£ A.t = 0, VI < k < n R , for some Ax <G AS 



if Hr is such 



Equivalently, for d min {HA,HB) = 0, the vectors hh,l < 

< n R should belong to the vector subspace [span(A;r)] 
for some Ax G AS riA+riB . In other words, the row space of 
the matrix Hr should be a subspace of the vector subspace 
[span(Aa;)] x , for some Ax G AS nA+nB , for d 2 min {H A ,H B ) 
to become zero. 

The channel fade coefficient matrix Hr is said to be a deep 
fade matrix, if d m i n (HA, H B ) = 0. The row space of the 
deep fade matrices are referred to as the deep fade spaces. 
The deep fade spaces are subspaces of the vector subspaces 
[span(Ax)] ± , Ax G AS nA+riB , which are referred to as the 
singular fade subspaces. Let J- denote the set of all singular 
fade subspaces, i.e., T = {[span(Aa;)] i : Ax G AS nA+7lB }. 

Note 1: Note that singular fade subspaces depend only on 
the set {Ax : Ax G AS nA+riB }, i.e., they are independent of 
nR. 

SISO two-way relaying as a special case: 

Consider the case when n A = ns = n R = 1. 
The row space of Hr, TZ(Hr) = k[HA H B ] T = 
k'[l jfjV i where k, k' G C. For this case, span(Ax) = 
c[Axa Axb) t , where c G C and Axa, Axb € AS. Hence, 
the singular fade subspaces [span(Aa:)] ± , Ax G AS 2 are of 
the form c'[l V \ where d G C. Note that for this 
case, TZ(Hr) and [span(Aa;)]" L are both one dimensional 
subspaces of the two dimensional vector space C 2 over C. 
Hence for TZ(Hr) to be a subspace of [span(A:r)] , it is 
necessary that TZ(Hr) = [span(Ax)] ± . It can be verified 
that U(Hr) = [span(Aa;)]- L if and only if -§a = =£r^- 
Hence, the effect of the MAI is totally captured by the 
ratio of the channel fade coefficients consistent with 

the results in (6), Ifl2l . Also, when the ratio referred to 

\ A 

as the fade state in lfF21 becomes equal to ~^ B A for some 
Axa,Axr G AS, the minimum distance of the effective 
constellation at the relay becomes zero. The complex numbers 
~: Axa , where Asa, Air G AS, were referred to as the 

Axb ' 5 ' 

singular fade states in lfl2l . 

From above, it is clear that the notion of singular fade 
subspaces subsumes in it as a special case the notion of 
singular fade states used for two-way relaying with single 
antenna at the nodes. 

Example 2: Consider the 2x2 system with BPSK signal 
set used at nodes A and B, with the bit mapped onto +1 
and the bit 1 mapped onto —1. The difference constellation of 



the BPSK signal set AS = {—2,0,2}. The vector subspace 
/ = [span([2 2 2 2} T )] 1 - = [span([-2 - 2 - 2 - 2} T )] 1 - is a 
singular fade subspace for this case. 

Broadcast (BC) phase: Let (xa,xb) G S nA+TlB denote the 
Maximum Likelihood (ML) estimate of (xa,xb) at R based 
on the received complex vector Yr, i.e., 



(xa,x b ) = arg min 

(x' A ,x' )eS«A+" 



\Y R - H a x'a - H B x\ 



B\- 



(3) 

Depending on Hr, R chooses a many-to-one map M Hr : 

S n A +n B where gl c C n a J s the signa j §et (of 

size between max {A/™- 4 , M nB } and M nA+nB ) used by R 
during the BC phase. Note that \S'\ should be at least 
max{M nA , M" B }, to transmit max{An,yi, Xn B } information 
bits. The elements in S nA+riB which are mapped on to the 
same complex vector in S' by the map J\A Hr are said to form 
a cluster. Let {£q, £2, Ct-i} denote the set of all such 
clusters. The formation of clusters for Hr is called clustering, 
and is denoted by C Hr . For simplicity, in the rest of the 
paper, the cluster Ck is denoted by the subscript k, where 
< k < t- 1. 

Example 3: For the 2x2 system with BPSK signal set, 



{[0 0] T , [0 1 1] T , [1 1 0] T , [1 1 1 1] T } 
[[0 1] T , [0 1 0] T , [1 1 1] T , [1 1 1 0] T ] 
> [0 10] T ,[0 111] T ,[1 0] T ,[110 1] T 
> [0 11] T ,[0 110] T ,[10 1] T ,[110 0] Ti 



(4) 



represents a clustering with four clusters. For two decoded 
pairs (xa,xr) and (x' A ,x' B ), the relay transmits the same 
vector from the signal set S' if [1^15] and [x' A x' B ] belong to 
the same cluster. For example, if R uses the clustering given 
in the vector transmitted during the BC phase will be the 
same, if the decoded pair during MA phase is ([0 0] T , [0 0] T ) 
or ([0 1] T , [0 1] T ), since [0 0] T and [0 10 1] T belong to 
the same cluster. 

The received signals at A and B during the BC phase are 
respectively given by, 



Y A = H'Xr + Z a ,Y b = H'Xr + Z t 



(5) 



where Xr = A4 Hr (xa,x b ) G S' is the complex vector 
transmitted by R. The fading matrices of size ua x ur and 
UbXtir corresponding to the R-A and R-B links are denoted 
by H' A and H' B respectively and the additive noises Za and 
Z B are CAf(0, a 2 ). 

In order to ensure that A (B) is able to decode B's (As) 
messages, the clustering C Hr should satisfy the exclusive law 
@, i.e., 



M Hr {x a ,x b ) ^ M Hr {x' a ,x b ),Vx a ^ x' A e S nA ,x B eS nB . 
M Hr {x a ,x b ) ^ M Hr (x a ,x' b ),Vx b ^x' B e S nB ,x A G S" A 



(6) 



Definition 1: The cluster distance between a pair of clusters 
Ci and Cj is the minimum among all the distances calculated 
between the points Haxa + H B x B and Hax' a + H B x' B G 



S R (H A ,H B ), where (x A ,x B ) G £s and (x' A ,x' B ) G 
The minimum cluster distance of the clustering C Hr is the 
minimum among all the cluster distances, i.e., 

d m in{C HR ) = min || H A (x A - x' A ) + H B (x B - x' B ) \\ . 

{x A ,x B ),{x' A ,x' B ')^s n A+ n B ■ 

M H R(x A ,x B )^M H R(x' A ,x' B ) 

The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during the 
BC phase is determined by the minimum distance of the signal 
set S 1 . Throughout, we restrict ourselves to optimizing the 
performance during the MA phase. For values of H B such that 
\hl Ax\ is small, VI < k < n R , for some Ax G AS nA+nB , 
dmin{Hn) is greatly reduced, a phenomenon referred as 
distance shortening. To avoid distance shortening, for every 
removable singular fade subspace, a clustering needs to be 
chosen such that the minimum cluster distance at every H B 
whose rows belong to that singular fade subspace is non-zero. 

For a singular fade subspace / G J 7 , let d m i n {C$ ,H B ) be 
defined as, 

d m in(C f ,H a ) = min \H A (x A - x' A ) + H B (x B - x' B ) |, 

{x A ,x B ),{x' A ,x' B )£S n A+"-B , 

(x A ,x B )^M f (x' A ,x' B ) 

where M.^ is the many-one map associated with the clustering 

cf. 

A clustering C$ is said to remove a singular fade subspace 
/ G J-, if the minimum cluster distance d m i n (C' , H B ) is 
greater than zero, for every H R such that 1Z(H R ) ^ /. If 
there are more than one clusterings which remove a singular 
fade subspace /, choose any one of them. 

It is important to note that certain singular fade subspaces 
cannot be removed. These are precisely the singular fade 
subspaces which are of the form [span(Ax)] ± , for which 

Ax = [0 nA Ax B ],Ax B G AS 7lB or 
Ax = [Ax A K B ], Ax A G AS n -\ 

and are referred to as the non-removable singular fade sub- 
spaces. The reason for this is as follows: The pair (x A ,x B ) 
and (x A ,x' B ) result in Ax = [0„ A Ax B \. But (x A ,x B ) and 
(x A , x' B ) cannot be placed in the same cluster since exclusive 
law given in (0 will be violated. Note that 0„ A +„ B is also a 
non-removable singular fade subspace, referred to as the trivial 
non-removable singular fade subspace. 

Note 2: The non-trivial non-removable singular fade sub- 
spaces for the SISO two-way relaying scenario are of 
the form [span([Aa; 0] T )} ± = span([0 1] T ) and 
[span([0 Ax] T )] ± = span([l 0] T ), where Ax e AS \ {0}. 
Irrespective of the signal set used at A and B, the non- 
trivial non-removable singular fade spaces are only two in 
number, since the vector subspaces [span([Ax 0] T )]' L are 
the same for all Ax £ AS \ {0}. In terms of the notion 
of singular fade states lfl2l . from the earlier discussion on 
SISO two-way relaying as a special case of MIMO two- 
way relaying, it follows that the above two singular fade 



subspaces correspond to the singular fade states infinity and 
zero respectively. Alternatively, the fact that these singular fade 
states are non-removable has been stated in [6] as follows: the 
distance shortening at H B /H A w (or H A /H B w 0) is 
inevitable. For SISO two-way relaying, the number of non- 
trivial non-removable singular fade subspaces remains two 
irrespective of the size of the signal set used at A and B. In 
contrast, it will be seen in Section II that for MIMO two-way 
relaying, the number of non-trivial non-removable singular 
fade subspaces increases with increasing size of the signal 
set used at A and B. 

Let Cjr = {C-f : / g J 7 } denote the set of all clusterings, 
which remove a removable singular fade subspace. 

Note 3: For every removable singular fade subspace, the set 
Cjr contains exactly one clustering which removes that singular 
fade subspace. The clusterings which belong to the set Cjr 
need not be distinct, since a single clustering can remove more 
than one singular fade subspace, as shown in the following 
example. 

Example 4: Consider the 2x2 system with BPSK signal 
set used at nodes A and B, with the bit mapped onto +1 
and the bit 1 mapped onto — 1. The difference constellation 
of the BPSK signal set is AS = {-2,0,2}. Consider 
the singular fade subspace, / = [span([2 2 2 2] T )] ± = 
[span([-2 - 2 - 2 - 2] T )] ± . The binary vectors 
[s Al s A2 s Bl s S2 ] T , [s' Al s' M s' Bi s' B2 } T G F| which result 
in the singular fade subspace / are [0 0] T , [1 1 1 1] T . 
The above two vectors need to be placed in the same cluster 
by the clustering which removes the singular fade subspace 
/. The clustering given in (|4]i in Example 3 removes the 
singular fade subspace /. Since [0 1 1] T and [1 1 0] T are 
in the same cluster, the clustering given in removes the 
singular fade subspace /' = [span([2 2 —2 — 2] T )] ± = 
[span([-2 - 2 2 2] T )]^ as well. 

In general, the row space of H B need not be a subspace 
of a singular fade subspace. In such a scenario, among all the 
clusterings which remove the singular fade subspaces, the one 
which maximizes the minimum cluster distance is chosen. In 
other words, for 7Z(H B ) ^ J 7 , the clustering C Hr is chosen to 
be C f , which satisfies d min (C f ,H R ) > d min (C f ' , H R ),Vf ^ 
/' G T . Since the clusterings which remove the singular fade 
subspaces are known to all the three nodes and are finite in 
number, the clustering used for a particular realization of the 
channel fade coefficients can be indicated by R to A and B 
using overhead bits. 

In J6|, a computer search algorithm called the Closest- 
Neighbour Clustering (CNC) algorithm was proposed for two- 
way relaying with single antenna at the nodes, which was 
used to obtain the network coding map that results in the 
best distance profile. The CNC algorithm was extended to 
the multiple antenna scenario in fTOl . The algorithm is run 
for a given H R . The total number of network coding maps 
which would result is known only after the algorithm is run for 
all possible realizations of H B which is uncountably infinite. 
Hence, the number of overhead bits required is not known 
beforehand. 



In contrast, the scheme proposed in this paper is based on 
the removal of singular fade subspaces. Since the number of 
singular fade subspaces is finite, the number of overhead bits 
required is upper bounded by the number of singular fade sub- 
spaces, which is known beforehand. In other words, the total 
number of network coding maps required is known exactly, 
which determines the number of overhead bits required. It is 
shown in Section III that the problem of obtaining clusterings 
which remove all the singular fade subspaces reduces to 
completing a finite number of partially filled Latin Rectangles, 
which totally avoids the problem of performing exhaustive 
search for an uncountably infinite number of values. 

The contributions and organization of this paper are as 
follows. 

• The structure and the exact number of non-removable and 
removable singular fade subspaces for M — PSK signal 
set (M any power of 2) are obtained analytically (Section 
II). It is shown that the fraction of the number of non- 
removable singular fade subspaces to the total number of 
singular fade subspaces tends to zero for large values of 
M (Section II). 

• It is shown that the requirement of satisfying the exclusive 
law is same as the clustering being represented by a 
partially filled Latin Rectangle (PFLR) and can be used to 
get the clustering which removes singular fade subspaces, 
by completing the PFLR (Section III A). 

• It is shown that the Latin Rectangles which remove the 
singular fade subspaces for the case when end nodes 
have unequal number of antennas, i.e., tia 7^ ns can 
be obtained from the Latin Squares which remove the 
singular fade subspaces for the case when the nodes have 
equal number of antennas n = m&x{riA, ns}(Section III 
B). 

• The singular fade subspaces which are removed by the 
conventional Exclusive-OR map are identified (Section 
III C). 

• It is shown that finding the network coding maps which 
remove all the singular fade subspaces reduces to finding 
a small set of maps. The entire set can be obtained from 
the small set by the notions of isotopic and transposed 
Latin Squares (Section III C). 

• The set of all Latin Squares which remove all the singular 
fade subspaces for the case when ua = «b = 2 and 
QPSK signal set is used at the end nodes is explicitly 
provided (Section IV). 

• It is shown that most of the Latin Squares which remove 
the singular fade subspaces for the n x n system, n > 
2, are obtainable from Latin Squares which remove the 
singular fade subspaces of the mxm system, where m < 
n (Section V). 

II. SINGULAR FADE SUBSPACES FOR 2 A -PSK 
SIGNAL SET 

In this section, the structure as well as the total number 
of singular fade subspaces is obtained for arbitrary 2 A -PSK 
signal sets. The points in the symmetric Af-PSK signal set 



are assumed to be of the form e jl - 2k+1)lT/M ,Q < k < M — 1 
and M is of the form 2 A , where A is a positive integer. 
For any M -PSK signal set, the set AS is of the form, 

AS = {0} U |2sin(7rn/M)e jfc27r / A/ |n oddj 

U {2 sin(7rn/M)e J '( fc27r / JU+ ' r '' M > |n even} , (7) 

where 1 < n < M/2 and < k < M - 1. 

In other words, the non-zero points in AS lie on M/2 
circles of radius 2 sm(irn/M ), 1 < n < M/2 with each circle 
containing M points. The phase angles of the M points on 
each circle is 2kir/M, if n is odd and 2kn/M + ir/M if n 
is even, where < k < M — 1. For example the difference 
constellation for QPSK signal set is shown in Fig. [3] 
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Fig. 3. Difference constellation for QPSK signal set 

The following Lemma is useful in finding the structure as 
well as number of singular fade subspaces. 

Lemma 1: For integers k\, k 2 , l\ and l 2 , where 

M 

1 < ki,k2,h,h < y,faf k 2 and h ^ l 2 , 

sin(fci7r/M) s'm(liir /M) 

— 5 (") 

sin(fc27r/M) sin(l2n/M) 

if and only if k\ = h and k 2 = h- 

Proof: See IH. ■ 
Recall from Section I that the singular fade subspaces are 
of the form [spanQAa;^ Ax^] T )] ± , where Ax A € AS nA 
and Ax B € AS" B . Let AxJ (Ax'g) denote the m th element 
of Axa (Axb)- Let 11,12, ■■■,iL be the ordered indices for 
which Ax]^ 7^ 0, 1 < k < L. Similarly, let Ji, J2, JL' 
be the ordered indices for which Ax J g 7^ 0, 1 < k < L' . Let 
'Pa ' 1 — — L anc ' 1 1 — ' — L' denote the phase angles of 



Ax l i and Ax^ respectively. Let the vector of length L+L'—l, 



■n) 



be referred to as the relative phase vector of Ax. 

Note that [span(Ax)]- L and [span(Aai')] can be the 
same for some Ax, Ax' G AS nA+TlB . A necessary and 
sufficient condition for [span(Ax)] ± = [span(Aa;')] ± , is 
that span(Ax) = span(Ax'), i.e., Ax = cAx', for some 
c G C. Equivalently, the conditions for [span(Ax)] 1 - = 
[span(Aa/)]- L can be stated as follows: 

• The location of the non-zero components is the same in 
the vectors Ax and Ax'. 

• The relative phase vector of Ax and Ax' are equal. 

• |Axi| = c|Ax' ; |, VI < i < n,A + tib, for some c G C. 
For 2 A PSK signal set, given the first condition, the third 
condition given above can be replaced by the condition given 
in the following lemma. 

Lemma 2: Let i\, £3, il be the ordered indices cor- 
responding to the non-zero components in Ax and Ax' (the 
location of the non-zero components is the same in the vectors 
Ax and Ax'). For 2 A -PSK signal set, \Axi\ = c|Ax-|, 
VI < i < ua + n>B, for some c G C, if and only if the 
magnitude of the non-zero components in Ax are equal and 
the magnitudes of the non-zero components in Ax' are equal, 
i.e., \Ax h \ = \Ax i2 \ = ...\Ax tL \ and |A<| = |A<| = 
...\Ax' iL \. 

Proof: When the condition given in the statement of the 
lemma is satisfied, clearly, |Axi| = c|Ax^|, VI < i < tia + 
n B . 

The proof of the "only if part is as follows. Let 1 < 
r\ < 7'2 < ... < tl < ua be the indices for which 
Ax^ fc ^ 0, Ax'\ k + 0, 1 < k < L. Similarly, let 1 < ji < 
32 < ■■■ < 3l < n B be the indices for which Ax^ ^ 
0,Ax' J b ^ 0, 1 < k < V . Since span(Ax) = span(Ax'), 
we have, 



A.j 



A.i 



Ax' 



Ax' 



(91 



From © it follows that the absolute value of the ratio of 
the points in the difference constellation AS are of the form 



,for some 1 < k, I < M/2. Equating the absolute 
values of the terms in (0 and from Lemma Q] it follows that 
the absolute values of the non-zero components in Ax need to 
be equal and the absolute values of the non-zero components 
in Ax' need to be equal. This completes the proof. 

■ 

Example 5: For the 2x2 system with BPSK signal set, for 
Ax = [0 2 2 - 2] T and Ax' = [0 - 2 - 2 2] T , the location of 
the non-zero components in Ax and Ax' are the same. The 
relative phase vectors of Ax and Ax' are equal. Also, we have, 

|Ax 2 | = |Ax 3 | = |Ax 4 | = 2,\Ax' 2 \ = \Ax' 3 \ = |Ax^| = 2. 
Hence, [span([0 2 2 - 2] T )] JL = [span([0 - 2 - 2 2] T )] ± . 



Example 6: For the 2x2 system with QPSK signal set, 
it can be seen from Fig. [3] that the non-zero points in the 
difference constellation lie on circles with radii \[2 and 2. 
For Ax = [0 V2 + jV2 s/2 - jV2 s/2 + jV2} T )} and 
Ax' = [0 V2 -jV2 V2] T )], the location of the non-zero 
components in Ax and Ax' are the same. The relative phase 
vectors of Ax and Ax' are equal. Also, we have, |Ax 2 | = 
|Ax 3 | = |Ax 4 | = 4,|Ax'| = |Ax'| = |Ax' 



the vector subspaces [span([0 
jVW)]^ and [span([0 y/2 



V2+3V2 V2- 

-3V2 VW)Y 



= 2. Hence 
-jV2 V2 + 
are the same. 



As mentioned in Section I, the singular fade subspaces 
can be classified in to removable and non-removable sin- 
gular fade subspaces. The non-removable singular fade sub- 
spaces are of the form [span([Ax J 4 0„ B ] T )] Axa G 
AS nA or [span([0„^ Ax B ] T )]^,Ax B G AS nB . The 
removable singular fade subspaces are of the form 
[span([Ax^ Ax^] T )] , Axa ^ 0„ A , Ax B 7^ 0„ B . 

The following lemma gives the total number of non- 
removable and removable singular fade subspaces. 

Lemma 3: For M -PSK signal set (M any power of 2), 

• the total number of non-removable singular fade sub- 
spaces is given by, 



E 



M 



M 
~2 

n B 



+ 1 



M 



M 



+ 1 



M 



the number of removable singular fade subspaces of the 
form [span(Ax)] ± with 2 < k < tia + n B non-zero 
components in Ax is given by, 



rn A + n B \ 
\ k J 



\ k ) 



(7) 



M 



M 
T 



M 



Hence the total number of removable singular fade 
subspaces is, 



n A +n I 

E 

fc=2 



where 



m A + n B ^ 
\ k ) 



(7) 



CD 



+ 1 



is defined to be zero if b > a. 



Proof: See Appendix A. ■ 

Example 7: Consider the MIMO two-way relaying system 
with ua = 2 and n B = 2. From Lemma [3] the non-removable 
singular fade subspaces for BPSK are 4 in number and are the 
ones numbered from 1 to 4 in Table I. For QPSK signal set, 
there are 28 non-removable singular fade subspaces. For BPSK 
and QPSK signal sets, the number of removable singular fade 
subspaces of the form [span(Ax)] ± , where Ax has k non-zero 
entries are as given below: 



No. 


Singular fade subspace 


1 


[span([0 2 2]' J ')] x =[span([0 -2 -2] T )] X 


2 


[span([0 2 -2]' J ')] x =[span([0 -2 2]' J ')] X 


3 


[span([2 2 0]' J ')] x =[span([-2 -2 0] T )] X 


4 


[span([2 -2 0]' J ')] x =[span([-2 2 0]' J ')] X 


5 


[span([0 2 2]' J )] x =[span([0 -2 -2]' J )] X 


6 


[span([0 2 -2]' J ')] x =[span([0 -2 2]' J ')] X 


7 


[span([2 2 0]' J ')] x =[span([-2 -2 0] T )] X 


8 


[span([2 -2 0] x )] x =[span([-2 2 0] x )] x 


9 


[span([2 2] x )] x =[span([-2 -2] x )] x 


10 


[span([2 -2]' J )] x =[span([-2 2]' J )] X 


11 


[span([0 2 2 0] T )] x =[span([0 -2 -2 0]' J ')] X 


12 


[span([0 2 -2 0]' J ')] x =[span([0 -2 2 0]' J )] X 


13 


[span([0 2 2 2]' J ')] x =[span([0 -2 -2 -2]' J ')] X 


14 


[span([0 2 2 -2]' J ')] x =[span([0 -2 -2 2] T )] X 


15 


[span([0 2 -2 -2]' J ')] x =[span([0 -2 2 2]' J ')] X 


16 


[span([0 2 -2 2] x )] x =[span([0 -2 2 -2] x )] x 


17 


[span([2 2 2] x )] x =[span([-2 -2 -2]' J )] x 


18 


[span([2 2 -2] T )] x =[span([-2 -2 2] T )] X 


19 


[span([2 -2 -2] r )] x =[span([-2 2 2]' J ')] X 


20 


[span([2 -2 2]' J ')] x =[span([-2 2 -2] T )] X 


21 


[span([2 2 2 F) x = spanf -2 -2 -2F)1 X 

L r VL J /J L r VL ^ 1/1 


22 


[span(2 2 -2]' J ')] x =[span([-2 -2 2] T )] X 


23 


[span([2 -2 -2]' J ')] x =[span([-2 2 2]' J ')] X 


24 


[span([2 -2 2]' J ')] x =[span([-2 2 -2] T )] X 


25 


[span([2 2 2 0] x )] x =[span([-2 -2 -2 0] J )] x 


26 


[span(2 2 -2 0]' J )] x =[span([-2 -2 2 0] x )] x 


27 


[span([2 -2 -2 0]' J )] x =[span([-2 2 2 0] T )] X 


28 


fspan([2 -2 2 F)F = spanf —2 2 -2 0F)F 

[ J r \ L J /J L r VL V J /J 


29 


[span([2 2 2 2F) x = span( -2 -2 -2 -2F)F 

L Jr VL J L Jr VL J /J 


30 


[span([2 2 2 -2 J ) x = span( -2 -2 -2 2l' J ')l x 

L I VL J /J L r VL J /J 


31 


[span([2 2 -2 -2 F) x = span( -2 -2 2 2F)1 X 

L r VL J /J L Jr VL J /J 


32 


[span([2 2 -2 2]' J )] x =[span([-2 -2 2 -2]' J ')] X 


33 


[span([2 -2 2 2]' J )] x =[span([-2 2 -2 -2]' J ')] X 


34 


[span([2 -2 2 -2]' J )] x =[span([-2 2 -2 2] r )] x 


35 


[span([2 -2 -2 2]' J )] x =[span([-2 2 2 -2] x )] x 


36 


[span([2 -2 -2 -2]' J )] x =[span([-2 2 2 2]' J )] X 



TABLE I 

SINGULAR FADE SUBSPACES FOR THE 2x2 SYSTEM WITH BPSK SIGNAL 
SET, WHERE THE SINGULAR FADE SUBSPACES NUMBERED 1 TO 4 ARE 
NON-REMOVABLE AND THE REST ARE REMOVABLE 



k 


No. of removable singular fade subspaces 


BPSK 


QPSK 


2 


8 


48 


3 


16 


448 


4 


8 


960 


Total 


32 


1456 



The removable singular fade subspaces for BPSK signal set, 
which are 32 in number, are the ones numbered from 5 to 36 
in Table I. 

From Lemma [3] it can be seen that the number non- 
removable singular fade subspaces is 0(M 2 max {™A,"B}-i^ 
while the number of removable singular fade subspaces is 
0(M 2 ( nA+nB > 1 ). Hence, the number of non-removable sin- 
gular fade subspaces is a small fraction of the total number 
of singular fade subspaces and the fraction tends to zero for 
increasing values of M. 



III. The Exclusive Law and Latin Rectangles 

For the two-way relaying scenario, with single antenna at 
the nodes, with signal sets of equal cardinality used at the 
end nodes, it was shown in [12] that all network coding maps 
satisfying the exclusive law are representable as Latin Squares. 
In this section, we establish the connection between Latin 
Rectangles and network coding maps satisfying the exclusive 
law, for the MIMO two-way relaying scenario. 

Definition 4: 01] A Latin Rectangle L of order Ni x N 2 
on the symbols from the set Z t = {0, 1, • • • ,t — 1} is an 
Ni x N2 array, in which each cell contains one symbol and 
each symbol occurs at most once in each row and column. A 
Latin Rectangle of order N x N is called a Latin Square of 
order N. 

Let the points in the A/ -point signal set used for transmis- 
sion at the nodes be indexed by the elements of the set TLm = 
{0, 1,2,... , M — 1}. Consider an M nA x M nB array at the 
relay with the rows (\columns) indexed by the n^-tuple^n^- 
tuple) [xA 1 ,XA a ,..- ,XA nA \ (\[xb 1 ,xb 2 ,--- ,x BnB }) denot- 
ing the complex vector transmitted by node A (\B). Our aim 
is to form clusters from the slots in the M nA x M nB array 
such that the exclusive law is satisfied. To do so, we will 
fill in the slots in the array with the elements of set Z t and 
the clusters are obtained by taking all the row-column pairs 
G Z^/ , j G Zjf such that the entry in the (i,j)-th 
slot is the same symbol from Z t for a cluster. The specific 
symbols from 7L t are not important, but it is the set of clusters 
that are important. Now, it is easy to see that if the exclusive 
law need to be satisfied, then the clustering should be such 
that an element in a row (column) cannot be repeated in the 
same row (column). Thus all the relay clusterings which satisfy 
the exclusive law form Latin Rectangles. Hence, we have the 
following: 

For MIMO two-way relaying, every relay clustering which 
satisfy the exclusive law forms a Latin Rectangle and vice 
verse. 

With this observation, the study of clustering which satisfies 
the exclusive law can be equivalently carried out as the study 
of Latin Rectangles with appropriate parameters. 

A. Removing singular fade subspaces, Singularity-removal 
Constraints and Constrained Latin Rectangles 

Consider a singular fade subspace / G F. Let (k, l)(k', I') £ 
l?lf x Z^f be the pairs which result in Ax such that 
[span(Ax)] - 1 = /. If (k, I) and (k',1 1 ) are not clustered 
together, the minimum cluster distance will be zero, for all Hr 
such that 1Z(Hr) < f. To avoid this, those pairs should be 
in the same cluster. This requirement is termed as singularity- 
removal constraint. So, we need to obtain Latin Rectangles 
which can remove singular fade subspaces and with minimum 
value for t. Therefore, initially we will fill the slots in the 
M nA x M" fl array such that for the slots corresponding to a 
singularity-removal constraint the same element will be used to 
fill slots. This removes that particular singular fade subspace. 
Such a partially filled Latin Rectangle is called a Constrained 
Partial Latin Rectangle (CPLR). After this, to make this a 
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Fig. 4. Constrained Partial Latin Square corresponding to the singular fade 
subspace [span([2 2-2 - 2] T )] 
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Fig. 5. Latin Square corresponding to the bit-wise XOR map for BPSK 
signal set 



Latin Rectangle, we will try to fill the other slots of the 
partially filled CPLR with minimum number of symbols from 
the set Z t . 

Example 8: Consider the 2x2 system with BPSK signal 
set used at the end nodes. Consider the singular fade subspace 
[span([2 2 - 2 - 2H] 1 - = [span([-2 - 2 2 2} T )}^. 
The singularity-removal constraint for this singular fade 
subspace is {([0 0], [1 1]),([1 1],[0 0])}. The constrained 
partial Latin Square for this case is shown in Fig. |4] 
The clustering which removes this singular fade subspace, 
given in in Example 3, can also be represented as a 
Latin Square shown in Fig. [5] For the singular fade subspace 
[span([0 — 2 — 2 2] T )] ± , the singularity removal constraints are 
{{([0 0] , [0 1]) , ([0 1] , [1 0])} , {([1 0] , [0 1]) , ([1 1] , [1 0])}} . 
The constrained partial Latin Square and the filled Latin 
Square which removes this singular fade subspace are shown 
respectively in Fig. [6] and Fig. [7] 

B. Obtaining Latin Rectangles from Latin Squares 

In this section, it is shown that the Latin Rectangles which 
remove the singular fade subspaces for the x ns system 
can be obtained from the Latin Squares, which remove the 
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Fig. 6. Constrained Partial Latin Square corresponding to the singular fade 
subspace [span([0 - 2-2 2] T )] ± . 
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Fig. 7. Latin Square that removes the singular fade subspace [span([0 
2-2 2] T )] X . 



singular fade subspaces for the n x n system, where n 
max{)u, kb}, by removing certain rows or columns. 



> riA, if [span([Ax A 
subspace for the tia x tib system, then 
1 is a singular fade subspace 



„TlT 



)]- 



is a sin- 



Ax A ' Ax' B ] r )]- 



For n B 
gular fade 

[ s Pan([0^ B _„ A 

for the tib x tib system. Similarly, for ua > Ub, if 
[spanQAa^ Ax^] 71 )] 1 - is a singular fade subspace for the 
Ua x tib system, then [span([AxJ 0n A -n B Ax B \ T )] 1 ~ is a 
singular fade subspace for the ua x ua system. 

For a Latin Square L of order n, let Lu- r A denote the Latin 
Rectangle of order r x n obtained by taking only the first r 
rows of L. Similarly, let £[ : ,i :c ] denote the Latin Rectangle 
obtained by taking only the first c columns of L. 

Lemma 4: For ns > n,A, if the Latin Square L removes the 
singular fade subspace [span([0^ B _ nA Ax^Ax^})] 1 - for the 

removes the 



-n A •''yl 

Us xtij system, the Latin Rectangle Li 1:M n A 



singular fade subspace [span([Ax^ Aa 



A l-lu, bS )] for the ha x n B 
system. Similarly, for tia > riB , if the Latin Square L removes 
the singular fade subspace [span([Ax^ 0„ A _ nB Ax^} 71 )} 1 - for 
the ua x ua system, the Latin Rectangle L^.a-.m^b] removes 
the singular fade subspace [span([Ax^ Ax^] 71 )]^ for the 
ua x tib system. 

Proof: Consider the case when ns > ua- Let / denote 
the singular fade subspace [span([0^ B _ AxJ Ax]^])] 1 - 
for the riB x ub system. Let Ax^ = [0„ B _„ A Ax A ] T . The 
singularity removal constraints for / are of the form, 



{([ X A[l:n B -n A ] X A[fi B -HA+l:riB]l > X b)j 



([ X A[l:n B -n A ] X ' B [n B -n A + l:n B ] ! X ' B ) } i 



Where XA[ nB -n A +l:n B ]] - x 'B[n B 

x' B = Ax B and Y.A[u nB - nA ] € S 

, the singularity removal constraints are 



n A + l:n B ] — Ax^, X# 

™.ForxAri:„, 



On^ — n 



{([0[i :n B -n A ] x A[TtB-ni+l:ns]l i x b)' 

([0[1:„ b -„a] X ' B[n B -n A + l:n B ]' X> b)}- ( 10 ) 

All the cells given in the constraints ( fTOb belong to the 
rows 1 to MriA- Since the Latin Square L which removes 
[span([0^ B _ IlA Ax^ AxU])] x satisfies the constraints in ( TTOl . 
the Latin Rectangle Lu-M n AA satisfies the constraints 

{([ X A[ nj j- ni + l: ni! ]] T > X B)i B[n B -n A +l:n B ]' X ' b)}> 

which are the singularity-removal constraints corresponding 
to the singular fade subspace [span([Ax^ Ax 1 ^)] 1 - for the 
ua'x ns system. ■ 

Example 9: Consider the 1x2 system with BPSK sig- 
nal set. The Latin Rectangle which removes the singular 
fade subspace [span([ —2 — 2 2] T )] ± shown in Fig. [8] 
is obtained by taking only the first two rows of the Latin 
Square in Fig. [7] which removes the singular fade subspace 
[span([0 - 2 - 2 2] T )] ± of the 2 x 2 system. 

From Lemma |4] it follows that the network coding maps 
for the case when tia ns can be obtained from network 
coding maps for the scenario in which the number of transmit 
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Fig. 8. Latin Rectangle that removes the singular fade subspace [span([— 2- 
2 2f)]-L. 



antennas at A and B n = max{n^, tib}- Hence, in the rest of 
the paper it is assumed that tia = Ub = n and the network 
coding maps to be obtained are Latin Squares. 

Note 4: If the maximum entry in the filled n x n Latin 
Square from which the Latin Rectangle is obtained is greater 
than n — 1, the obtained Latin Rectangle also can have entries 
greater than n— 1. Since the performance during the BC phase 
is dependent on the number of distinct entries in the filled Latin 
Rectangle, this adversely impacts the performance during the 
BC phase. In that case, obtaining the Latin Rectangle by direct 
completion of CPLR, instead of obtaining it from the Latin 
Square, may lead to a better performance. 

C. Some Special Constructions of Latin Squares 

Recall that the rows and columns of the Latin Squares are 
indexed by vectors which belong to Zjy-. By bit- wise XOR 
of two such vectors, it is meant the vector obtained by taking 
the bit-wise XOR of the individual components of the two 
vectors, after decimal to binary conversion. Every cell in the 
Latin Square corresponding to the bit-wise XOR mapping is 
filled with the decimal equivalent of the bit-wise XOR of the 
row index and the column index. 

Consider the singular fade subspaces 
[span([Ax y iAa;B] T )]" L , which satisfy the condition that 
AxAi = ± Axs t , VI < i < n. Let T± denote the set of such 
singular fade subspaces. In the following lemma, it is shown 
that bit-wise XOR mapping removes all the singular fade 
subspaces which belong to T±. 

Lemma 5: When the user nodes use 2 A -PSK constellations, 
the singular fade subspaces which belong to the set T± are 
removed by bit-wise XOR mapping, for all A. 

Proof: By definition, for / = span^A^^A^] 1 ^) G 
J 7 ±, the non-zero locations of Ax a and Axb should match. 
Let L be the number of non-zero components in Axa 
and Axb- Without loss of generality assume, Axj = 
[Ax Jl Axj 2 ... Axj L O n - L ] T , J e {A, B}, (if the non- 
zero components of Axa and Axb appear in any other 
order, the indexing given for the transmit antennas can be 
permuted to get the assumed ordering). By definition, any 
/ G T± should satisfy, Axa { — ±Aa;_B»,l < i < L. Let 
XAt = ki, 1 < i < L and XBi = h t l < i < L. For 



{[fci k 2 ...k L h l 2 ...l L ] T ,Wi k' 2 ...k' L l[l' 2 



■1'l\ T 



} 



to be a singularity-removal constraint, since AxAi = ± Axb ; , 
it follows that, 



e M = ± e M 



,l<i<L. (11) 
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Fig. 9. Latin Square that removes the singular fade subspace 
[span([-2 2 -2] T )] X . 



Equating the magnitude and phase terms in (fTTT i results in 
the following possibilities: 



ki + fcj — li + ij + li, 

Ajj k^ — lj or/uj 



k'i — h h + 



Solving ( fT2l and ( TT3l , we get ki 



li , ki 



(12) 

M. (13) 

= V, or ki = 



/' 



M 7. _ 1,1 i M 
2 , H — ft-i -r 2 



Similarly for the case when Axa 
should satisfy are, ki 

,k' 



/' /■ 



, the conditions 

k': or ki = L + 



ki,li, k[ and 
M 

2 

XOR map places in the same cluster those vector pairs (k, I) 



l[ + 4p . Hence, it needs to be shown that the bit- wise 



and (fc',Z'), for which ki = li,k[ 
k[ + 4f, or ki = I'^k = k^ or k z 

M 



M 
2 ' 1 
i M 



l[ or ki = l[ + 

ii + f,K = l 'i + f> 

l'i + f,yi<i<L. 



ki — l'i j li — k^ or ki — li H Y^'i 

The Latin Square corresponding to the bit-wise XOR map 
has the following properties: 

• ki + ki = 0. Hence, it places vectors with fe, = h, k[ = l[ 
in the same cluster. 

• ki + k = ki — 4f + li — 4f . Hence, it places vectors with 
ki = l'i + 4p , li = k'i + 4r m the same cluster. 

• h + h = li + ki. Hence, it places vectors with ki = 
I'ifli = k[ in the same cluster. 

M _ ; i / M 



ki + ki — #■ = li + k — M- . Hence, it places vectors with 
i + 4p , k'i = l'i + M- in the same cluster. 



k 



Consider the case when Axa- = Axb 



This completes the proof. ■ 
Example 10: Consider the 2x2 system with BPSK signal 
set. The singular fade subspaces 6-9 and 30-37 (12 in total) 
given in Table 1 belong to F± and are removed by the bit-wise 
XOR map. 

Definition 2: A Latin Square L T is said to be the Transpose 
of a Latin Square L, if L T (i,j) = L(j,i) for all i,j G 
{0,1,2,.,M-1}. 

Lemma 6: If the Latin Square L removes the singular fade 
subspace [span([Ax^ Aa;^)] 7 "]- 1 , then the Latin Square L T 
removes the singular fade subspace [span([Ax^ Ax^)] T ] ± . 

Proof: The singular fade subspace 
[spanQAa;^ Aa;^)] 1 "]^ can be viewed as the singular 
fade subspace [span([Ax^ Ax^)] 71 ]- 1 with the users A 
and B interchanged. Interchanging the users is the same 
as interchanging the row and column indices, i.e., taking 
transpose. ■ 

Example 11: For a 2 x 2 system with BPSK signal set, 
since the Latin Square given in Fig. [7] removes the singular 
fade subspace [span([0 —2 —2 2] T )] ± , from Lemma [6] its 
transpose shown in Fig. |9]removes the singular fade subspace 
[span([-2 2 -2] T )] ± . 
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Fig. 10. Latin Square obtained from XOR map by the permutation of row 
indices 
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Fig. 11. Latin Square which is same as the one shown in Fig. 1101 
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Fig. 12. Latin Square that removes the singular fade subspace [span([0 2 — 
2 2H]X. 
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Fig. 13. Latin Square that is same as the Latin Square in Fig. 1121 



Definition 3: |[TT1 Two Latin Squares L and L ' (using the 
same symbol set) are isotopic if there is a triple (f,g,h), where 
/ is a row permutation, g is a column permutation and h is 
a symbol permutation, such that applying these permutations 
on L gives L'. 

Consider a vector Ax = [Ax^ Ax|J] t , where Ax a and 
Axb are obtained by the applying the permutations a a and 
Ob on Ax a and Axb respectively. Equivalently this can 
be viewed as applying the permutations a a and ffg on the 
indices of the transmitting antennas at tia and n b respectively. 
Since the rows and columns of the Latin Squares which 
remove the singular fade subspaces are indexed by the vectors 
transmitted by nodes A and B respectively, applying the 
permutations a a and as on the components of the row and 
column indices of the Latin Square which removes the singular 
fade subspace [span(Ax)] x results in an isotopic Latin Square 
which removes the singular fade subspace [span(Ax)] x . This 
is stated as the following lemma. 

Lemma 7: If a Latin Square L removes the singular fade 
subspace [span(Ax)] ± , the Latin Square V obtained by apply- 
ing the permutation a a on the components of the row indices 
and the permutation ub on the components of the column 
indices of L removes the singular fade subspace [span(Ax)] J -. 

Example 12: Consider the 2x2 system with BPSK signal 
set. As seen in Example 10, the XOR map given in Fig. [5] 
removes the singular fade subspaces 6-9 given in Table 1. 
Permuting the components of the row indices of the Latin 
Square in Fig. (i.e., the first component becomes the second 
component and vice verse), results in the Latin Square Fig. 
[TOl which is the same as the Latin Square shown in Fig. QT| 
From Lemma [7] it follows that the Latin Square in Fig. [TTI 
removes the singular fade subspaces 9-13 given in Table 1. 

Consider the case when A/-PSK signal set (M any power 
of 2) is used at the end nodes. Consider two singular fade 
subspaces [span(Ax)] 1 - and [span(Ax)] 1 - which are such that 
the absolute values of the components of Ax and Ax are 
equal, i.e., they differ only in the relative phase vectors. Let 
[4>A 4>b\ T and \4>' T a </>'b] T be the relative phase vectors of Ax 
and Ax respectively, where <f>A and (j)' A are of length n — 1, 
and <\>b and <fi' B are of length n. Let Acf>A = 4>a ~ 4>'a an d 



A(j) B = 4>b-<I>'b- Also let A^Ai = -^r and A<j} Bi = -3^- 
Lemma 8: For 2 A -PSK signal set, let L denote the 
Latin Square which removes the singular fade subspace 
[span(Ax)]- L . The Latin Square L 1 which removes the singular 
fade subspace [span(Ax)] ± can be obtained from L as follows: 
To the i th component of all the row indices of L add 
^AjjVl < i < n modulo M and to the i th component of 
all the column indices of L add fcg,;,Vl < i < n modulo M, 
to obtain the Latin Square L' . 

Proof: Rotating the phase of the i th component of Ax^ 
by an angle A £ v , can be viewed equivalently as rotating the 
phase of the signal set used by the antenna rii at node A, by the 
same angle. In the Latin Square L that removes [span(Ax)]" 1 -, 
this change can be effected by adding fc^ modulo M to the 
i th component of the row indices. By a similar reasoning, ks, 
needs to be added to the column indices of L to obtain V . 
This completes the proof. ■ 
Example 13: Consider the 2x2 system with BPSK signal 
set. As discussed in Example 8, the Latin Square shown in Fig. 
[7]removes the singular fade subspace [span([0— 2— 22] T )] ± = 
[span([0 2 2— 2] T )]~ L . From Lemma [8] adding 1 modulo 2 to 
the second component of the row index in the Latin Square in 
Fig. [7] we get the Latin Square shown in Fig. [12] which is the 
same as the one in Fig. Qj] which removes the singular fade 
subspace [span([0 2 -2 2] T )] ± . 

The usefulness of Lemmas |6]|8] is that the set of all Latin 
Squares which remove all the singular fade subspaces can be 
obtained from a small set of Latin Squares. This is illustrated 
for the 2x2 system in the following subsection. 

D. Sufficient number of Latin Squares to be obtained for 2x2 
system with M-PSK signal set 

Consider the 2x2 system with A/-PSK (M any power 
of 2) signal set used at A and B. For this case, the 
singular fade subspaces are of the form [span(Ax)] x where 
Ax = [Ax A t Axa 2 Ax Bl Axb 2 ] t e AS 4 . Let k be the 
number of non-zero components in Ax. The following lemma 
gives the sufficient number of Latin Squares from which the 
Latin Squares which remove all the removable singular fade 



subspaces can be obtained. 

Lemma 9: For the 2x2 system, with 2 A -PSK signal set, 
the sufficient number of Latin Squares from which the Latin 
Squares which remove all the removable singular fade sub- 
spaces of the form [span(Ax)] ± , with k non-zero components 
in Ax can be obtained are given by, 



k 


No. of Latin Squares 


2 


M* M , i 
,-fl , 4 + 1 


3 


AI , M 2 M , i 
Ifi "T" 8 2 ~r 1 


4 


A/ 4 | M i 11M' Z 3M | 1 
128 32 32 8 "T ^ 



Proof: Case 1: k=4 
For this case, from Lemma [3] there are 
(if) — TT sm g u l ar fade subspaces. 

For vectors Aa; and Ax such that |Axj| = |Axj|,Vl < i < 
n i.e., Ax and Ax differ only in their relative phase vectors, 
from Lemma|7] the Latin Square which removes [span(Ax)]- 1 - 
can be obtained by row and column permutations. This reduces 
the total number of Latin Squares to be obtained by a factor 
M 3 . In other words the total number of Latin Squares to be 



M 
2 



1 



{{Axi,Ax 2 } G AS 2 



M(M+2) 
8 

such 



choices 
that the 



obtained is upper bounded by 
Consider the set AS 2 - 

2 ) + ( 1 ) 

for the pair {Axi,Ax2}. If Ax and Ax 
unordered pairs {AxaAxb} and {Ax^Axs} are the same, 
from Lemma [8] the Latin Square which removes the singular 
fade subspace [span(Ax)] 1 " can be obtained from the one 
which removes [span(Ax)]- 1 . 

Hence the number of Latin Squares to be obtained is equal 
to the number of ways of choosing two elements from the set 
AS? minus the number of cases which results in the same 
singular fade subspaces. The total number of such possibilities 
is. 



A/(A/+2) 
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Af(A/+2) 
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128 
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Case 2: k=3 For this case, from Lemma [3] the number of 



;f) 3 



M 
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1 



M 1 



removable singular fade subspaces are 
Considering only those singular fade subspaces for which | Ax\ 
is distinct, the total number of Latin Squares to be obtained 
for this case is less than (4m — 4p + 1 . 

The number of Latin Squares to be obtained is equal to the 
number of ways of choosing an element from the set AS? Q 
and an element from the set AS \ {0} minus the number of 
cases which results in the same singular fade subspaces. The 
total number of such possibilities is, 



M(M + 2) 
8^ 



Case 3: k = 2 
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For this case, from Lemma [3] there are [-^| 4p + l]M 3 

singular fade subspaces. The total number of Latin Squares to 
be obtained is equal to the total number of ways of choosing 
two elements from A5\{0} minus the number of cases which 
result in the same singular fade subspaces. The total number 
of such possibilities is, 



M 
T 



1 = 



8 



M 
T 



1. 



Note 5: Lemma [9] provides only the sufficient number of 
Latin Squares from which all the Latin Squares can be 
obtained. The actual number can be lesser than the number 
given in Lemma [9] since the same Latin Square can remove 
more than one singular fade subspace. For example, from 
Lemma it follows that the bit-wise XOR map removes more 
than one singular fade subspace. 

IV. ILLUSTRATIONS FOR THE 2x2 SYSTEM WITH 
QPSK SIGNAL SET 

Consider the 2x2 system with QPSK signal set. The non- 
zero points of the difference constellation of the QPSK signal 
set AS lie on two circle with radii s/2 and 2. The four points 
on the circles with radius y2 have phase angles anc [ 
the 4 points on the circle with radius 2 have phase angles 
(2m M 1 ^ ,0<^<3. 

As seen in Example 7, there are 1456 removable singular 
fade subspaces, out of which 960 have k = 4. From Lemma 
[9] it follows that the set of Latin Squares which remove all 
these 960 singular fade subspaces can be obtained from 5 Latin 
Squares, which remove the following singular fade subspaces: 

h = [span([\/2 + jV2 V2 + 3V2 V2 + 3V2 V2 + jV^})}^, 

/a = [span([\/2 + jy/2 V2 V2 + jV2 V2 + jV^])]" 1 , 

h = [span([\/2 + jsfl V2 V2 V2])] x , 

U = [span([V2 + j\/2 V2 V2 + j\/2 V2])]- 1 and 

fs = [span([V2 + jV2 V2 + jV2 V2 v^])]^. 

From Lemma [5] the singular fade subspaces f\ and fi are 
removed by the bitwise-XOR map, given in Fig. |17(a)| given 
in Appendix B. The Latin squares which remove the other 
three singular fade subspaces and /s are given in Fig. 

[T7(b)lfT7(d)l in Appendix B. 

For k = 3, the Latin Squares which remove all the 448 sin- 
gular fade subspaces can be obtained from the Latin Squares 
given in Fig. |17(e)fT7(i)| in Appendix B, which remove the 
following 5 singular fade subspaces: 

fe = [span([0 y/2 + 3V2 V2 + jV2 V2 + j\/2] T )] ± , 

f 7 = [span([0 V2 + jV2\/2\/2 + jV2] T )] i , 

fs = [span([0 V2 + 3V2 V2 V2] T )] ± , 

fa = [span([0 V2 + jV2 V2 + jV2 \/2] T )] ± and 

f w = [span([0 V2V2 + jy/2 \/2] T )] ± . 
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(a) The Latin Squares L (b) The Latin Square L ® L 

and Z 

Fig. 14. Example illustrating the notion of Cartesian product of Latin Squares 



For k = 2, there are 48 removable singular fade sub- 
spaces. From Lemma [9] the Latin Squares which re- 
move all these 48 singular fade subspaces can be ob- 
tained from the 2 Latin Squares which remove the singular 
fade subspaces fn = [span([0 y/2 \/2] T )] i - and /12 = 
[span([0 V2 + jV20 V2] T )] ± . From Lemma EJ it follows that 
the XOR map given in Fig. |17(a)| in Appendix B removes fn. 
The Latin Square that removes /12 is given in Fig. 17(j)| in 
Appendix B. To sum up, the Latin Squares which remove all 
the 1496 removable singular fade subspaces, are obtainable 



from the 10 Latin Squares given in Fig. [I7]in Appendix B. 

V. OBTAINING LATIN SQUARES FOR THE n x n 
SYSTEM FROM LATIN SQUARES OF LOWER 
ORDER 

In this section, it is shown that most of the Latin Squares 
which remove the singular fade subspaces of the n x n system, 
n > 2 are obtainable from the Latin Squares which remove 
the singular fade subspaces of m x m system, where m < n. 

Definition 4: For two vectors y and z of length 2a and 
26 respectively, the compound vector of y and z, denoted as 
comp(y, z) is the vector of length 2a + 2& given by, 

r T T T T iT 

[V[l:a] Z [l:b] V[a+l:2a\ Z [b+l:2b}\ ■ 

For example, for two vectors y = [1 2 3 4] T and z = [0 5] T , 
comp(j/, z) = [1 2 3 4 5] T . 

For a Latin Square L, let Lu.j^.j] denote the (j — i + 1) x 
(I — k + 1) array obtained by taking only the i th to j th rows 
and k th to I th columns of L. Let L + c denote the Latin Square 
obtained by adding integer c to all the cells of L. Let m&x{L} 
denote the maximum among all the integers filled in the cells 
of the Latin Square L. 

Definition 5: The Cartesian product of the two Latin 
Squares L of order M a and L of order M , denoted as (L<S)L), 
is the Latin Square of order M a+b for which 

(i®i)[(i_i)Af"+l:iAf»,0'-l)M«+l:3Af] = L + L(i,j) (max{L} + 1) , 

where 1 < i,j < M b . 



For example, the Cartesian product of the Latin Squares L 
and L of order 2 1 and 2 2 respectively shown in Fig. |14(a)| is 
the Latin Square of order 2 3 shown in Fig. |14(b)| 

Consider two vectors Ax <G AS 2a and Ax £ AS 2b . The 
vector subspaces [span(Ax)] ± and [span(Ax)] ± are singular 
fade subspaces for the a x a and b x b systems respec- 
tively. The following lemma shows that the Latin Square 
which removes all the singular fade subspaces of the form 
[span(comp(Ax, kAx))]^, k € C can be obtained by taking 
the Cartesian product of the Latin Squares which remove 
[span(Ax)] ± and [span(Ax)]^. 

Note 6: For x G AS 2b and k e C such that kAx G S 2b ', 
the vector subspaces of C 2b [span(A.T)] ± and [span(fcAx)]- 1 - 
are the same. But the vector subspaces of C 2a+2b 
[span(comp(Ax, AS))]- 1 and [span(comp(Aa;, kAx))}- 1 , k ^ 
1 are different. 

Lemma 10: Let L and L respectively denote the Latin 
Squares of order M a x M a and M b x M b , which remove 
the singular fade subspaces [span(Ax)] ± of the ax a system, 
and [span(Ax)] ± of the b x b system. The Latin Square 
L ® L removes all the singular fade subspaces of the form 
[span(comp(Ax, kAx))]^ , k £ C. 
Proof: For the ax a system, let 



{{[x lA X2 A ...X aA ] , [x lB X 2a ■■■Xa B \) , 

A A X 2 A ■■■ x 'a A \> \ X 1 B x 2j 



'-])} 



denote a singularity removal constraint for the singular fade 
subspace [span(Aa;)] ± . Similarly, for the b x b system, let 

{([yi A V2 A ■■■Vb A ] , [yi B V2 A ■■■Vb B \) > 

( [y'i A V2 A --- y' bA ] , [y'i B y r 2 A • • • y' bB ] ) } 

denote a singularity removal constraint for the singular fade 
subspace [span(Ax)] ± . Since span(Ax) = fcspan(Ax), 

{ ( A k yiA ■ ■ • k Vb A ] > k y? A ■ ■ ■ k VbB ] ) ' 

{[ky[ A ky' 2A ...ky' bA ],[ky' lB ky' 2A . . . kyQ)} 

is also a singularity removal constraint for the singular fade 
subspace [span(Ax)]- L . 

It can be verified that the Cartesian product of the 
Latin Squares which remove the singular fade subspaces 
[span(Ax)] ± and [span(Ai)] ± satisfy the Cartesian product 
of their constraints. In other words, the Latin Square (L Cg> L) 
satisfies the constraint, 

{{[xi A x 2a ...Xa A ky lA ky-2 A ...ky bA \ , 

[xi B x 2a ...Xag ky lB ky 2A ...ky bB \) , 
( [x' lA x' 2a . . . x' aA ky[ A ky' 2A . . . ky' bA ] , 

[x' lB x' 2a ...x' aB ky' lB ky' 2A . . . ky' bB \) } , 

which is a singularity removal constraint for the singular fade 
subspace [span(comp(Aa;, kAx))} 1 - of the (a + b) x (a + b) 
system. This completes the proof. ■ 



From Lemma [TUJ it follows that Latin Squares for those 
singular fade subspaces for the n x n system, which are 
expressible as [span(comp(Aa;, kAx))] , where [span(Aa;)] 
and [span(Ax)] ± are removable singular fade subspaces for 
a x a and 6x6 systems respectively, for some choices of 
a < n and 6 < n such that a + b = n, are obtainable from the 
Latin Squares for the a x a and 6x6 systems. 

Example 14: Consider the 3x3 system with BPSK signal 
set. Consider the singular fade subspace [span([2 — 2 2 2 2 — 
2] )]-*-. The Latin Square L shown in Fig. |14(a)| removes the 
singular fade subspace [span([2 2] T )^ of the lxl system and 
the Latin Square L shown in Fig. 1 14(a)] removes the singular 
fade subspace [span([— 222 — 2])} 1 - of the 2 x 2 system. Since 
[2 -2222-2]= comp([2 2], [-2 2 2 - 2]), the Latin Square 
L<g)L shown in Fig. 1 14(b)] removes the singular fade subspace 
[span([2 - 2 2 2 2 - 2])] x of the 3 x 3 system. 

From Lemma[3]it follows that the number of non-removable 
singular fade subspaces for the n x n system with A/-PSK 
signal set is of order Af 4n_1 and the singular fade subspaces 
of the form [span(Aa-)] x , where Ax has all the components to 
be non-zero contributes to the maximum number of singular 
fade states which is of order M 4n ~ l . It is easy to verify that a 
vector Ax with all the components to be non-zero can always 
be written of the form comp(Axi, Aa^), where Axi and Ax2 
are of lengths a and 6, where a,b < n and a + b = n, such that 
[span(A.Ti)] ± and [span(Ax2)] ± are removable singular fade 
subspaces for the ax a and 6x6 systems respectively. Hence, 
it follows that most of the Latin Squares which remove the 
singular fade subspaces for the n x n system are obtainable 
from Latin Squares of lower order. 

VI. SIMULATION RESULTS 

The proposed scheme is based on the removal of all the 
singular fade subspaces, i.e., a minimum cluster distance 
greater than zero is ensured for all realizations of the channel 
fade coefficient matrices. Also, it is ensured that the number 
of clusters in the clustering, which is the same as the size 
of the signal set used during the BC phase is minimized. 
As a reference scheme, we consider the case when bit-wise 
XOR network code is used at R, irrespective of the channel 
conditions. Since XOR network does not remove all the 
singular fade subspaces, the proposed scheme is expected 
to perform better than the pure XOR network code based 
scheme, which is confirmed by the simulation results. All the 
simulation results presented are for the case when QPSK signal 
set is used at the end nodes. The number of antennas at all 
the three nodes are taken to be two. 

The additive noises at the nodes are taken to be of unit 
variance. It is assumed that signal sets of equal energies are 
used at all the nodes, which is defined to be the Signal to Noise 
Ratio (SNR). Fig. [T5l shows the SNR vs BER performance for 
the proposed scheme and the for the case when XOR network 
code is used irrespective of the channel condition, for the case 
when each one of the entries of the matrices Ha, Hb, H' a 
and H' B follow i.i.d. Rayleigh distribution, with unit variance. 
Fig. Q~6] shows a similar plot for the case when each one of 
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Fig. 15. SNR vs BER for different schemes for QPSK signal set for a 
Ravleiah fading scenario with channel variances dB 
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Fig. 16. SNR vs BER for different schemes at node B for QPSK signal set 
for a Rician fading scenario with a Rician factor of 10 dB 

the entries of the matrices Ha, Hb, H' a and H' B follow i.i.d. 
Rician distribution with a Rician FactoiQ of 10 dB and unit 
variances. From Fig. [15] and Fig. [16] it can be seen that the 
diversity order is two, for the proposed scheme as well as 
the scheme which uses the conventional XOR network code. 
Also, it can be seen that for the Rayleigh and Rician fading 
scenarios considered, at a BER of 10~ 4 , the proposed scheme 
provides a gain of 1.5 dB and 2.5 dB respectively over the 
conventional XOR network code based scheme. 

VII. Discussion 

An adaptive network coding scheme using Latin Rectangles 
for the wireless two-way relaying scenario was proposed. The 
scheme was based on the removal of a finite number of vector 
subspaces referred to as the removable singular fade subspaces 
by proper choice of network coding maps. For a particular 

'Rician factor is the power ratio between the line of sight and scattered 
components. 



channel realization, among all the network coding maps which 
remove the singular fade subspaces, the one which maximizes 
the minimum clustering distance is chosen. In this way, the 
set of all possible channel realizations, which is C nA+nB for 
the tia x riB system, is quantized in to a finite number of 
regions, based on which one of the network coding maps 
obtained optimizes the performance. Such a quantization was 
obtained analytically in 1131 for the lxl system. Obtaining 
the quantization for a general tia x ?i_b system is an interesting 
problem for further investigation. 
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Appendix A 
PROOF OF LEMMA[3] 

Number of non- removable singular fade subspaces: 

Consider the case when the non-removable singular fade 
subspaces are of the form [span([Aa;yi Ong] 71 )]^. Let k 
be the number of non-zero components of Ax a- Assume 
that the relative phase vector of Axa is fixed. Since the 



points in AS lie on 4j- circles, there are (4=f) possibil- 
ities for the absolute values of the non-zero components 
of Axa- But out of these, if the absolute values of all 
the components of Ax a, are equal, from Lemma [2 the 
resulting singular fade subspaces are the same. Subtracting 
all such cases which are 4^ m number and adding one 



to account for all such cases, results in (4f) — 4f + 1 
singular fade subspaces. From Lemma [2] given the same 
absolute values of the non-zero components of Ax a, distinct 
singular fade subspaces results for distinct relative phase 
vectors of Axa- Since there are M k ^ 1 distinct possibilities for 
the relative phase vector, there are 



singular fade subspaces for the case considered 



with k 

non-zero components. Summing over all possible values 
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k=l \ k I 
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sin- 



of k, we have 
gular fade subspaces for the case considered. Similarly, 
for the case when the non-removable singular fade sub- 
spaces are of the form [span([0„_ 4 Axb] 71 )]^, there are 



3 (?) 



1 = 1 



M 
2 



1 



M 1 singular fade subspaces. 



Number of removable singular fade subspaces: 

Let k be the number of non-zero entries in Ax = 
[Ax a Ax^] T . Since Axa ^ 0„ A and Axb ^ 0„ B , k should 
be at least 2. Without loss of generality, assume that tib > tia- 
Case 1: 2 < k < ua 

The k non-zero components of Ax can be chosen in 
( nA ~l nB ) ways. But this also includes the cases when all 
the k non-zero components occur in Axa and Axb = 0„ B 
or vice verse. Subtracting out those cases, the k non-zero 
components of Ax can be chosen in (™ 4 +™ B ) - (™ fe 4 ) - (™ fe B ) 
ways. For each one of those possibilities, following an ap- 
proach similar to the one used in the proof of Lemma 



k-l 



[3] it can be shown that there are [(4/ ) fe - 4f + l] M 
singular fade subspaces. Hence, for this case there are 

[m - (Y) - (T)] [(¥)* - f + 1] M"- 1 possible 
singular fade subspaces in total. 

Case 2: tia < k <ub 

The k non-zero components of Ax can be 
chosen in (™ A fc™ B ) — (™jf) ways, since the case 
when all the k non-zero components occur in 
AxB,i.e., Axa = 0„ 4 , needs to be excluded. Hence, there 

m [(«a+-b) _ (naj] _ M + j_ j M k-1 possibilities 

in total for this case. 

Case 3: tib < k < ua 

The k non-zero components of Ax can 

be chosen in (™' 4 fc™ B ) ways. This results in 

( nA t nB ) [(¥) fc - TT + !] sin g ular fad e subspaces in 

total for this case. 

For 2 < k < tia + tib, summing up the number of 
possibilities obtained in the above three cases and defining 
m to be zero for b > a, the number of removable singular 
fade subspaces is as given in the statement of the lemma. This 
completes the proof. 



Appendix B 

LATIN SQUARES FOR THE 2x2 SYSTEM WITH 
4-PSK SIGNAL SET 

For the 2x2 system with 4-PSK signal set, the Latin Squares 
from which the set of Latin Squares which remove all the 
removable singular fade subspaces are obtainable, are given 
in the next three pages. 
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(a) Latin Square that removes the singu 


ar tade subspaces fi, /4 and /n 
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(b) Latin Square that removes the singul 


ar tade subspace /a 
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(c) Latin Square that removes the singular fade subspace fs 
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Fig. 17. Latin Squares that remove different singular fade subspaces 
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Fig. 17. (Contd.) Latin Squares that remove different singular fade subspaces 
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Fig. 17. (Contd.) Latin Squares that remove different singular fade subspaces 



